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We discuss selectron pair production in e~e~ scattering. These processes can 
only occur via t-channel neutralino exchange, provided that the neutralinos are 
. Majorana fermions. We concentrate on the processes e~e~ — > e^e^, grgr — ► 

e ~Xi e ~X2 ~* e ~Xi e_ Xi//> hi which a complete determination of the final state 
momenta is possible without using the selectron masses as input. The exper- 
imental reconstruction of the selectron masses in this decay channel provides 
clear evidence of the Maiorana character of the neutralinos, which is confirmed 

cn ■ 

by the analysis of the electron energy spectrum. 
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1 Introduction 

43 : 



Supersymmetry (SUSY) is one of the most interesting and natural extensions of the 
standard model of fundamental interactions (SM). In addition to providing a solution 
to many of the SM problems, supersymmetric theories allow a natural connection to 
high energy scales, and thus to more fundamental theories of particle physics [T |l2|l3] . If 
SUSY is present at the weak scale, it is expected that it will be discovered at LHC jlj, 
or even at the second run of Tevatron [H]. After the new particles have been detected, 
it must be verified that they are indeed the superpartners of the SM fields. It will be 
crucial to measure their quantum numbers and, for the case of neutral gauginos, con- 
firm that they are Majorana fermions. Moreover, masses, mixings, couplings and CP 
violating phases must be measured in a model-independent way, so that the Lagrangian 
parameters can be determined and the SUSY relations for the couplings verified. These 
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tasks are not easy to accomplish, because in general it will be difficult to separate the 
different sectors, and many processes giving the same experimental signatures will be 
simultaneously present. In order to disentangle the different processes and obtain pre- 
cise measurements, a high energy e + e~ collider like TESLA is essential to complement 
the LHC capabilities [BJ. In fact, the TESLA design offers many advantages for SUSY 
studies, with the possibility of beam polarisation and the option of e~e~, and 77 
scattering. 

Here we focus on the determination of the Majorana nature of the neutralinos at 
TESLA. In previous works [7,8.9], it has been pointed out that the neutralino character 
can be tested in the processes e + e~ — > XxX^ — * with % = 2,3,4 and I = 

e, //. These processes give an experimental signature of two leptons of opposite charge 
plus missing energy and momentum 1 , and the energy distributions of the final state 
charged leptons are sensitive to the Dirac or Majorana nature of the decaying neutralino 
Here we follow a different approach, studying selectron pair production in e~e~ 
collisions [HUH]. In contrast with e + e~ scattering, the processes e - e~ ee (with 
e = ei,e.R) are only mediated by diagrams with t-channel neutralino exchange. The 
Majorana nature of the neutralinos is essential for the nonvanishing of the transition 
amplitudes, as can be clearly seen in Fig. [TJ We note that in some SUSY models it 
is possible that one or more pairs of Majorana neutralinos combine to form Dirac (or 
pseudo-Dirac) fermions T3|. This happens when there are two Majorana neutralinos 
with the same mass (or approximately the same mass, for the pseudo-Dirac case) 
and opposite CP parities. In this situation, the contributions of these two Majorana 
neutralinos to selectron pair production are (almost) equal in modulus with opposite 
signs, and cancel (or nearly cancel) in the amplitude. 



A.? 



\ / 
\ / 
\ / 
/ \ 
/ \ 
/ \ 



Figure 1: Feynman diagrams mediating the processes e e ->ee. 



1 We assume that the first neutralino x\ is the lightest supersymmetric particle (LSP), which is 
stable if R parity is conserved. Since the \\ is neutral, colourless, and weakly interacting, it escapes 
detection, producing a signature of missing energy and momentum. 
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The dominant decay mode of the selectrons is e — > e~xi, leading to a final state 
with two electrons plus missing energy and momentum. The observation of this signal 
at sizeable rates is a clear indication of selectron pair production, although in principle 
it could also originate from other new physics processes. In this decay channel, the 
masses of the selectrons cannot be directly reconstructed, due to the presence of two 
unobservable neutralinos in the final state. However, the selectron masses can be 
precisely determined, for instance with a measurement of the total cross section at 
energies near threshold (TOIITT] . In the continuum, the study of the electron energy 
spectrum shows that the produced particles are scalars, and the selectron masses can 
be extracted from the end points of this distribution jTHj. Furthermore, it is possible to 
construct the minimum kinematically allowed selectron mass |14| . whose distribution 
peaks at the actual selectron masses and hence yields another measurement of these 
quantities. The coincidence of all these measurements provides evidence for selectron 
pair production, and disfavours the interpretation of the e~e~ signal as originated by 
other new physics process. 

With the same purpose, in this paper we consider a different decay channel which 
offers the advantage that the 4-momenta of all final state particles can be determined, 
and thus the selectron masses can be directly reconstructed event by event. This 
can be achieved restricting ourselves to tifih and e^eR production, and selecting the 
channel where one of the selectrons decays to e~\\ an d the other selectron decays via 
e — » e~%2 — ¥ e ~Xiff> with // a or qq pair. Although the latter is a relatively 

rare decay mode, if the selectrons are light enough to be produced at TESLA, the high 
luminosity available will allow us to obtain an observable signal in most cases. Using 
this channel alone, the direct reconstruction of the masses and the analysis of the two 
electron energy spectra (one resulting from e — > e~Xi and the other from e — > e~x%) 
provides clear evidence for selectron pair production. Since the reconstruction method 
introduced here allows the determination of the momenta of all final state particles, it 
can also be used for further studies of the angular distributions in the production and 
decay of the selectrons. In particular, it provides an independent direct verification 
that the selectrons are scalar particles. This analysis will be presented elsewhere [T5] . 

This paper is organised as follows: In Section El we discuss the production and 
subsequent decay of selectron pairs in e~e~ collisions, and how these processes are gen- 
erated. In Section El we briefly describe the model, outlining some important features, 
and we select two illustrative scenarios for and qq final states. In Section 0] we 

describe how the reconstruction of the selectron masses is achieved and we present our 
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results. Section [HI is devoted to our conclusions. In Appendix El we collect the mass 
matrices and Lagrangian terms which are relevant for our work. The benefits of beam 
polarisation for these processes are displayed in Appendix iBl 

2 Production and decay of selectron pairs in e~e~ 
scattering 

The reconstruction of the selectron masses is only possible for the production of two 
identical selectrons, e~e~ — > eifit,, crCr (this will be justified in the detailed discussion 
of Section IH). For these processes we choose the decay channels with one selectron 
decaying via e — > e^Xi an d the other e — > e~x\ ~> e Xiff- I n addition, we consider 
mixed selectron production e~e~ — > enei, with tR — > e~Xi, &l — ► e~x\ ~^ e ~Xiffi 
since this is the main background to the two processes of interest 2 . Each of the processes 

e~e 
e~e 
e~e 

is mediated by 56 Feynman diagrams, depicted in Figs. El and El Although for clarity 
the diagrams for the decay of X2 are shown separately, in our computations all the 
diagrams are summed coherently. We only consider final states with // = fi + fi~,qq, 
and in the case of qq we sum uu, dd, ss, cc and bb production, without flavour tagging. 
Concerning the remaining channels, in X2 ~ * Xi e+e ~ the multiplicity of electrons in 
the final state makes it difficult to identify the electron resulting from the decay of the 
X*2- I n X2 ~> Xi u v t ne presence of four undetected particles in the final state yields 
too many unmeasured momenta to allow their kinematical determination. The same 
happens in X2 Xi T+T ~ 1 because each of the r leptons decays producing one or two 
neutrinos that escape detection. 

Compared with tRth production, the other SUSY and SM backgrounds are less 
important. For instance, the process e~e~ — > eR&Rff, with both selectrons decaying 
to e~Xii nas a smaller cross section, since it has the same number of vertices as eR&L 
(both amplitudes are proportional to g 6 ) and does not have the enhancement of the 
X2 near its mass shell. The main SM background is e~e~ — ► e~e~vi>ff, including 

2 In the same process e - e~ — ► BrCl, the channel Sl — > e~Xi, en — > e^\2 ~ > e ~Xiff nas a much 
smaller cross section because of the smaller branching fractions (see Section . 



e L e L -> e~x±e~X2 ~* ^Xi^XiJJ 

e R e R -> e~Xie~X°2 ^ e~x?e~x?/7 

e R e L -> e"Y?e~Y2 -> e~Y?e~Y?// (1) 
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Figure 2: Feynman diagrams for the processes in Eqs. (IT}. The 4 neutralinos are ex- 
changed in the t channel, and the shaded circles stand for the 7 sub-diagrams mediating 
the decay of separately shown in Fig. El 




(a) (b) 




(c) (d) 

Figure 3: Feynman diagrams for the decay of mediated by Z bosons (a), neutral 
Higgs bosons (b), and left-and right-handed scalar fermions (c and d). 

resonant e~e~~ — > W~e~vff and e~e~ — > e~e~Zff production. This process has a 
small cross section and is separable from the e^ex, &r&r signals, which have larger 
values for the missing energy and momentum. One might also consider the process 
e~e - — > W~W~ ff mediated by heavy Majorana neutrinos. The cross section for the 
latter is very constrained by present limits on neutrinoless double beta decay [IE], and 
this process could also be separated from selectron pair production, using the missing 
energy and momentum. 
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We calculate the full matrix elements for the resonant processes in Eqs. (|T]) with 
HELAS [17J, using the Feynman rules for Majorana fermions given in Ref. ^B]- The 
inclusion of next-to-leading order corrections is not necessary in this work, since it 
only increases slightly the cross sections [Hi]. We assume a centre of mass (CM) energy 
Ecm = 2-E = 500 GeV for the e~e~ collisions, and an integrated luminosity of 100 fb _1 . 
In order to take into account the effect of initial state radiation (ISR) on the effective 
beam energies, we convolute the differential cross section evaluated for fractions xi, x 2 
of the beam energies with "structure functions" Disr(^i) and -Disr(^2), namely 

-i 



da 



da{x 1 E, x 2 E)Dis R {x 1 )Di SR (x 2 )dx 1 dx2 . 



(2) 



The explicit expression for -Disr used is [21]] 
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7 is the Euler constant, a = 1/137 the fine structure constant, s = E^ M and m e the 
electron mass. The effect of beamstrahlung is taken into account in a similar way using 
the "structure function" [2T] 
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with z = (y/3) 3 / 4 . We take N = 0.555, re = 0.0693 [21] for our evaluations. 

In order to simulate the calorimeter and tracking resolution, we perform a Gaussian 
smearing of the energies of electrons (e), jets (j) and muons (//) using the specifications 
in the TESLA Technical Design Report |23j 
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where the energies are in GeV and the two terms are added in quadrature. We apply 
"detector" cuts on transverse momenta, p? > 10 GeV, and pseudorapidities \rj\ < 2.5, 
the latter corresponding to polar angles 10° < 6 < 170°. We also reject events in which 
the leptons and/or jets are not isolated, requiring that the distances in (rj, 6) space 
AR = a/ Ar/ 2 + AO 2 satisfy AR > 0.4. We do not require specific trigger conditions, 
and we assume that the presence of charged leptons with high transverse momentum 
will suffice. Finally, for the Monte Carlo integration in 8-body phase space we use 
RAMBO [24J. 

3 Overview of the model and selection of scenarios 

The aim of the present work is to demostrate that the selectron masses can be suc- 
cessfully reconstructed with the method here proposed. It is nevertheless important to 
show that this can be done in realistic scenarios. In this section, we outline the main 
features of the model, and investigate the behaviour of the selectron production cross 
sections and decay rates as a function of the free parameters of the model. 

In our analysis we consider the minimal supersymmetric standard model (MSSM), 
where a minimal number of fields is introduced and i?-parity is conserved. For sim- 
plicity, we assume an underlying minimal supergravity (mSUGRA) [22) framework. At 
the grand unification scale, Mqut — 2 x 10 16 GeV, relations of universality among 
the several soft breaking terms allow the model to be described by five independent 
parameters: the common soft scalar mass (m ), the unified gaugino mass (m 1( / 2 ), the 
universal trilinear coupling (A Q ), the ratio of the two neutral Higgs vacuum expectation 
values (tan/3 = Vijv\) and the sign of the bilinear Higgs term, sign(/x). In our work 
we do not address supersymmetric CP violation, thus we assume that mi/ 2 and A 
are real. The low-energy Lagrangian is derived by means of the renormalisation group 
equations, which are numerically solved to the two-loop order using SPheno (2H]. In 
addition, SPheno is used to obtain the masses and mixings of SUSY particles, as well as 
some decay widths. The mass matrices for sfermions and neutralinos and the required 
interaction terms are presented in Appendix El 

It is important to clarify here some issues, which are relevant for the subsequent 
discussion of selectron production and decay. The first concerns flavour and chirality 
mixing in the slepton sector. For the case of the first two generations, the LR term 
in the slepton mass matrix (Eq. (|A.3j) ) is significantly suppressed by the electron and 
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muon Yukawa couplings. The smallness of LR mixing for selectrons and smuons is 
fairly insensitive to the particular value of A , and for simplicity we choose A = in 
our numerical analysis. This typically translates into 9g ~ 0(1CT 4 ) and 6^ ~ 0(1O~ 2 ) 
(cf. Eq. (|A.6jl ). so that clPr and [ilI^r mixing can indeed be neglected. Therefore, 
for the first two generations the mass eigenstates are the chiral (e, (e, states. 
Moreover, leptonic interactions are assumed to be flavour-conserving, so only selectrons 
are produced in e~e~ collisions. 

As afore mentioned, we assume that the lightest neutralino is the LSP. Depending on 
the underlying model, the Xi can be gaugino- or higgsino-like, or a balanced admixture 
of the latter interaction eigenstates. Here we shall only consider cases where x\ is 
gaugino-like and, more specifically, dominated by the bino component, a situation 
common to a vast region of the mSUGRA parameter space. In our case, the second 
neutralino, i s dominated by the wino component, and thus its coupling to the 
right-handed selectron e# is very small. This fact has a significant impact in our 
analysis, since in this situation the e# mainly decays to e~ x®, even when the e~ x\ 
channel is kinematically allowed, so that the crCr process in Eq. JH) has a very small 
cross section. We emphasise that the reconstruction method proposed in this paper is 
not sensitive to the nature of the lightest neutralinos. However, the cross sections of 
the processes studied will be significantly affected, given that the relevant x'J-fermion- 
sfermion couplings are suppressed if the lightest neutralinos are higgsino-like. 

The numerical studies presented in this section do not aim to be a complete scan 
of the available parameter space, but rather an illustrative analysis of the impact of 
the variation of the mSUGRA parameters on the selectron production and decay cross 
sections. In the following, we will simultaneously try to satisfy the direct and indirect 
bounds on the masses of the particles [27J and still obtain an adequate spectrum for 
sleptons and neutralinos, so that the processes in Eq. JT} can be observed at TESLA 
with a CM energy of 500 GeV. 

We begin the selection of the scenarios by considering the dependence of the studied 
signals on m . This is shown in Fig. 01 for e^ei,, e-R&L and e^eR production (full, 
dashed and dot-dashed lines) and for the fi + fi~ and qq channels (black and grey lines, 
respectively). We have taken mu2 = 220 GeV, tan/3 = 10, and /i > 0. The shaded area 
on the left corresponds to values of m excluded by the current experimental bounds 
on mg R [2Z| • I n addition, the end points of the curves marked with a circle correspond 
to the values of m below which the lightest neutralino is not the LSP. In this section 
we consider unpolarised beams; the effect of longitudinal polarisation is discussed in 
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Figure 4: Cross sections (in fb) for e e — > ee — ► e X2 > e Xi e Xi// (with 
unpolarised beams) as a function of m , for m 1 / 2 = 220 GeV and tan/3 = 10. 



From this plot we notice that only for m > 140 GeV the e^ is heavy enough to decay 
into the second neutralino. The shape of the curves reflects the strong dependence on 
mo of the production cross sections and branching ratios for ex — > z~X% an d X2 — > Xiff- 
In the latter case, there are several competing Feynman diagrams [28,29], mediated by 
the Z boson and sfermions, that contribute to the decay of the X2 (the contribution 
of the diagrams with a Higgs boson is small for this value of tan/3). Three different 
regimes can be distinguished: 

• For m < 120 GeV, the second neutralino decays predominantly to charged lep- 
tons, X2 ~^ Xil + l~- The decay amplitudes are dominated by diagrams (c) and 
(d) in Fig. El with the exchange of on-shell right-handed sleptons, in particular 
the decay %\ — > i s dominated by fiR exchange. The contribution of dia- 
gram (a) with an on-shell Z boson is less important, due to the smallness of the 
ZX2X1 coupling. In this region of the parameter space, the decays to XiQQ are 
very suppressed, not only because of the small coupling of the Z boson to Xi and 
X21 but also due to the heavy squark masses, nig > 400 GeV. 

• For m ~ 130 GeV, the e# and jiR are heavier than the x°, while the tr is lighter. 
Then, the X2 decays almost exclusively to Xi T+T ~ ■, mediated by an on-shell tr, 
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diagrams (c) and (d) in Fig. 01 (As remarked in Section [2J in the r + r~ channel 
the selectron masses cannot be reconstructed, and thus we do not consider these 
final states in the plots presented throughout this section.) 



• For m > 140 GeV, the right-handed sleptons (including f R ) are heavier than 
X°, and the Z-exchange diagram in Fig. El dominates, yielding a large branching 
ratio for x° — * XiQQ- 

For m > 200 GeV, the production of e^L and eifiR is not possible with a CM 
energy of 500 GeV, and only tRtR pairs are produced. From inspection of Fig. 0] we 
select two values m = 80 GeV and m = 160 GeV, approximately in the centre of the 
first and third regions, that will serve to illustrate the reconstruction of the selectron 
masses for fi + fi~ and qq final states. The sets of parameters for these two mSUGRA 
scenarios are summarised in Table QJ Both of them are in agreement with the bounds 
arising from the b — ► 57 decay and they are in the region of the (mi/ 2 ,m ) plane 
favoured by — at the 2 a level [30]. The resulting selectron and neutralino masses 
and widths, as well as the relevant branching ratios, are collected in Table 01 for each of 
these scenarios. In scenario 1 the tR is lighter than the and the decay e# — > e~x\ 
is not possible. Hence, in this scenario we only consider the eifih signal. 



Parameter Scenario 1 Scenario 2 
m 1/2 220 220 
m 80 160 
A 
tan/3 10 10 
sign fi + + 



Table 1: Input parameters for the two mSUGRA scenarios to be considered in Section 
El The values of muz, m o and A are in GeV. 



For completeness, we investigate the dependence of the cross sections on the re- 
maining parameters. For scenario 1, and leaving all other parameters as in Tabled we 
plot in Fig. [HI the cross sections for different values of m^. The shaded area on the left 
corresponds to values of mi/2 that yield a Xi mass excluded by present limits |22j. Note 
that the depression occurring around m\/2 = 160 GeV has the same origin as the one 
observed in Fig. SI and corresponds to the dominance of the r + r~ decay channel. The 
decrease of the cross sections for large mut is caused by the reduction in the branching 



10 



Scenario 1 Scenario 2 

m gL 181.0 227.4 

T e - L 0.25 0.85 

m- eR 123.0 185.0 

T iR 0.17 0.58 

m,o 84.0 84.3 

Xi 

m,o 155.8 156.4 

T^o 0.023 1.50 x 10~ 5 

m,o 309.4 310.0 

I> 1.48 1.43 

m^o 330.4 331.1 

T^o 2.30 2.01 

Br(g L -> e~xi) 41.8 % 18.3 % 

Br(g L -> e-Xz) 20.6 % 30.8 % 

Br(g fl -> e~x°i) 100 % 99.7 % 

Br(£ R -> e"X2) - 0.3 % 

Br(x° -> 10.3 % 3.9 % 

Br(x§ -> ~ 69.2 % 



Table 2: Some relevant quantities in the two mSUGRA scenarios considered in Section 
EJ The masses and widths are in GeV. 

ratio for — > e~X2 and i s a l so a consequence of the smaller phase space available for 
selectron pair production. In fact, with m = 80 GeV, the production of e^e-h pairs is 
only kinematically possible when < 340 GeV. 

For scenario 2, we fix m = 160 GeV, tan/5 = 10, A = (as in Table HI), and 
plot the cross sections for different values of mu2 (see Fig. IH)). The shaded area on 
the left of this figure corresponds to values of m\/i for which m^o is excluded [2*7) . 
The sudden decrease of the cross sections for mi/ 2 > 250 GeV is caused by both the 
dominance of the t + t~ channel and the small phase space available for selectron pair 
production. For m = 160 GeV, the production of e^L pairs is kinematically forbidden 
with mi/2 > 270 GeV. The production of e^eR pairs is still allowed, but for m-i/2 > 260 
GeV the \2 i s heavier than the e.R and the decay £r — ► e~\2 is n °t possible. 

We also examine the dependence of the cross sections on tan (3 for scenarios 1 and 
2. Fixing m 1/2 = 220 GeV, m = 80, 160 GeV, A = (as in TableHJ, we plot in Figs.EI 
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Figure 5: Cross sections (in fb) for e e — > ee — > e Xi e X2 ~ ¥ e Xi e Xiff (with 
unpolarised beams) as a function of mi/2, for m = 80 GeV and tan/3 = 10. 




m 1/2 (GeV) 

Figure 6: Cross sections (in fb) for e~e~ — > ee — > e^Xi e xH ~ > e Xi e Xiff (with 
unpolarised beams) as a function of m^, for m = 160 GeV and tan/3 = 10. 



and |H1 respectively, the cross sections for different values of tan (5. In both figures the 
shaded areas on the left are excluded by the LEP bounds on the MSSM lightest Higgs 
mass 3 , and in Fig. the end points of the curves marked with a circle correspond 
3 Since theoretical calculations of m^o have an estimated error of 2-3 GeV [HI], we assume the 
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to values of tan/3 for which the LSP is no longer Xi- In both cases the regions of 
parameter space for which the process e~e~ — > eifiL —> e ~ Xi e ~Xi ff is observable 
range up to tan/3 ~ 22. 




tan P 



Figure 7: Cross sections (in fb) for e e — > ee — > e x? e X2 > e Xi e Xiff (with 
unpolarised beams) as a function of tan /3, for mi/ 2 = 220 GeV and m = 80 GeV. 

The study of Figs. 0HH1 reveals that for most regions of mSUGRA parameter space 
where tifih production is kinematically accesible, either the ji^ or the qq decay 
channels have observable cross sections. The only exceptions are narrow intervals in 
mo and mi/2, and also large tan (3 values, where the r + r~ channel completely dominates 
the X2 decays. On the other hand, the cross section for e~e~ — > crEr — > e~Xi e ~X2 ~^ 
e ~Xi e ~Xiff is always small, as a consequence of the small branching ratio for e# — > 
e~X2 (i n addition, this decay mode is not possible when m^ R < m^o)- Finally, we 
remark that, as shown in Figs. for larger values of m\n and/or m selectron 
pair production and, in general, the production of supersymmetric particles, is not 
kinematically allowed at the assumed CM energy of 500 GeV, and a higher energy 
collider (or the TESLA upgrade) would be required for studies involving a heavier 
SUSY spectrum, 
conservative bound m h a > 111 GeV. 
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tan P 

Figure 8: Cross sections (in fb) for e~e~ — > ee — ► e _ Xi e ~X2 — * e ~Xi e ~Xiff (with 
unpolarised beams) as a function of tan/?, for min = 220 GeV and m = 160 GeV. 

4 Reconstruction of the selectron masses 

To reconstruct the selectron masses we use as input the 4-momenta of the detected 
particles (the two electrons and the // pair), the CM energy and the Xi an d X° masses, 
which we assume known from other experiments [6lll3| . Remarkably, the values of both 
selectron masses are not needed anywhere. Let us first fix our notation by labelling 
the two electrons in the final state as "electron a" and "electron 6" (either choice is 
equivalent) and their momenta as p a = (E a ,p a ) and p b = (E b ,p b ). The selectrons 
to which they correspond are then called "selectron a" and "selectron 6", and the Xi 
resulting from the decay of the latter are also labelled as "a" and "i>", with momenta 
q a = (E' a ,q a ) and q b = (E' b ,qb), respectively. The momentum of the // pair is denoted 
as p f j=(E f j,p ff ). 

In general, it is necessary to have as many kinematical relations as unknown vari- 
ables in order to determine the momenta of the undetected particles. In our case, there 
are 8 unknowns (the 4 components of q a and q b ) and 8 constraints. From energy and 
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momentum conservation and the fact that the two Xi are on shell, we have 

E a + E b = E C m — (E a + E b + Eff) , 

Qa + Qb = i , 
2 2 

q a = m-o, 

ql = m 2 -o , (9) 

with $ the missing momentum. Another equation is obtained from the decay of the 
intermediate Since we do not know a priori which selectron decays to we m ust 
consider both possibilities. If we assume that it is the selectron "a", we have 

(la+Pff) 2 = m 2 -o [case (a)] , (10) 

or, assuming it is the selectron "6", 

(Qb+Pff) 2 = m 2 -o [case (b)] . (11) 

In order to set the last kinematical constraint, we make the hypothesis that in the e~e~ 
scattering two particles of equal mass are produced 4 . This implies that in one case we 
have 

E a + E' a + E ff - = E b + E' b [case (a)] (12) 

while in the other 

E a + E' a = E b + E' b + E ff [case (b)] . (13) 

It is worthwhile remarking that ISR, beamstrahlung, particle width effects and 
detector resolution degrade the determination of the \i momenta. The ISR and beam- 
strahlung modify the effective beam energies, and therefore the first two of Eqs. (JHI 
do not hold exactly. The finite detector resolution also implies that Eqs. (j9T [T3|) are 
only approximate. Moreover, the off-shellness of the selectrons and of the xl have 
a non-negligible influence on Eqs. (fTnHT3|l . As a result of all these effects, in some 
cases the solutions of these equations are imaginary, even when the x° i s "correctly" 
assigned to a selectron. In order to find a real solution in any circumstance, we force 
the discriminants of the second degree equations involved to be non-negative. 

Each set of equations - Eqs. (I9fl()ll2|) and Eqs. (|9I11I13|) - yields two solutions for 

q a and q b . Among these four solutions, we have to select one, which in turn gives a 

value for the reconstructed mass. To do so, we observe the following criteria: 

4 This applies for e^L and es&R production, but not for £r(Il, in which case the selectron masses 
cannot be reconstructed. 
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(I) We first eliminate the unphysical solutions yielding negative reconstructed selec- 
tron masses m| ec . 

(II) We then discard the solutions in which the reconstructed mass of the x%i ?n™o, is 
too different from the real value 5 : for the channel we require |ml e c — m^o| < 

5 GeV, and for the qq channel |m~ e c — m^ol < 10 GeV. 

(III) Among the remaining solutions, we take the one giving the smallest value of mf c . 
This choice provides the "correct" solution in most cases. If no solution is left 
after (II), the event is discarded. 

The above procedure determines the momenta of the two unobserved xli identifying 
which selectron has decayed to x° as well, allowing to reconstruct the selectron masses. 
The reconstruction also allows to distinguish between the electrons resulting from e — > 
e~xS) an d e — > e~X2- Let us define E\ and E 2 as the energies of the electrons produced 
in the decays e — > e~x1, and e — * e"X2) respectively. In the selectron rest frame, the 
electron energy is fixed by the kinematics of the 2-body decay, and since the selectrons 
are spinless particles, the decay is isotropic. Then, for e^e^, &r&r production the 
energy distribution of the electrons in the CM frame is flat, with end points at [T41I32] 



where /3 = y 1 — 4mJ £ J s. For mixed selectron production the electron energy spectra 

are flat as well, but the expressions of the end points are more involved. The results 
for the reconstruction are illustrated for scenarios 1 and 2. 

4.1 Scenario 1 

As seen from Fig. in scenario 1 the reconstruction of the mass is quite effective, 
and a sharp peak at the true value = 181 GeV (as given in Table is obtained. 
On the other hand, the distribution of the tRth events is rather flat. We notice that the 
criteria (I, II) in the reconstruction procedure operate as an efficient kinematical cut, 

5 Although by construction mi e c = m f o, forcing the system of equations to have a real solution 
relaxes this equality in some cases (see the preceeding paragraph). 




(14) 
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suppressing the unwanted e^e-i, production, while for eifih the effect is much smaller 
(cf. Table 03). Assuming an integrated luminosity of 100 fb _1 , the number of events 
collected in the detector is sufficiently large to allow the observation of a positive signal, 
with a mass distribution concentrated at mg L = 181 GeV. The use of negative beam 
polarisation P\ = P 2 = —0.8 improves these results, increasing the cl^l cross section 
by a factor of 3.24 and reducing the crCl cross section by a factor of 0.36 (see Table El 
and Appendix lB|. 

1.25 



1.00 



> 0.75 
u 
O 
tS 

°0.50 
0.25 

'loo 150 200 250 300 

m. (GeV) 

Figure 9: Reconstructed selectron mass in scenario 1, for unpolarised beams. 

Poo P- 

before after before after 

e L e L 3.25 2.88 10.52 9.32 

e R e L 1.66 0.47 0.60 0.17 

Table 3: Cross sections (in fb) for selectron pair production in scenario 1. We quote 
results before and after reconstruction, for unpolarised beams (-Poo) an d for Pi = P 2 = 
-0.8 (P__). 

The particles observed in the kinematical distribution in Fig. El are scalars, as can 
be easily deduced for instance from the angular distributions for the production and 
the decay [11]. The same can also be concluded from the analysis of the electron energy 
spectra, which also provides additional independent determinations of their mass. From 
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Scenario 1 
m m = 220 GeV 
m„ = 80 GeV 
10 



tan |3 = 




Eqs. ijTijl. the expected end points of the distributions are E™ m = 30 GeV, Ef 1 ^ = 166 
GeV and Ef n = 10 GeV, £ 2 max = 55 GeV. The kinematical distributions of E 1 and 
E 2 are shown in Figs. [TQl and [TTJ respectively. We have included the e^L signal, as 
well as the crCl background. The energy spectrum of the electrons is smeared by 
ISR, beamstrahlung, detector resolution and particle width effects. However, the end 
points can be clearly observed in both plots, further supporting the evidence for e^L 
production. 



0.08 



0.06 



0.04 



0.02 



Scenario 1 
m [/2 = 220 GeV 
m = 80 GeV 
tan P = 10 




100 
E, (GeV) 



200 



Figure 10: Kinematical distribution of E\ in scenario 1, for unpolarised beams. 



4.2 Scenario 2 

In this scenario, the clCl events produce a peak at mg L = 227 GeV, and the crCr 
events yield a tiny peak around vn^ R = 185 GeV, as depicted in Fig. [12J Noticeably, 
the distribution of the crCl events concentrates around 215 GeV. This behaviour is a 
result of the smaller ratio (mg L — rrie R ) / {rrie L + fng R ): in this scenario, the hypothesis 
of two particles produced with equal mass, used for the reconstruction, becomes more 
accurate. For the same reason, in this case this background is less reduced by the 
criteria (I, II) than in scenario 1. We observe that the e^L peak is broader than in 
scenario 1, due to the larger width and the less precise measurement of the energy 
for jets. Additionally, the distribution is slightly concentrated on smaller rrie R values, 
as a consequence of the procedure used for the reconstruction. 
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Figure 11: Kinematical distribution of E 2 in scenario 1, for unpolarised beams. 
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Figure 12: Reconstructed selectron masses in scenario 2, for unpolarised beams. 

The use of negative beam polarisation enhances the e-i&L signal and reduces the 
crCl background as in scenario 1, practically eliminating e^e^ production (see Table EJ). 
Conversely, the use of positive beam polarisation P\ = P 2 = 0.8 enhances crCr by a 
factor of 3.24, suppressing eh&L and tRtL by factors of 0.04 and 0.36, respectively, and 
making the 6r6r peak clearly visible in this shown in Appendix [Bj 
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Poo 




P_ 


- 


P+ 


+ 




before 


after 


before 


after 


before 


after 




6.52 


5.71 


21.14 


18.50 


0.26 


0.23 


e R e R 


0.41 


0.34 


0.02 


0.01 


1.34 


1.10 


eRe L 


5.41 


2.63 


1.95 


0.95 


1.95 


0.95 



Table 4: Cross sections (in fb) for selectron pair production in scenario 2. We quote 
results before and after reconstruction, for unpolarised beams (Poo), f° r Pi — P2 — —0.8 
(P__) and for P 1 = P 2 = 0.8 (P++). 

The examination of the energy spectra in Figs. PCTl and 1141 confirms the evidence of 
production. From Eqs. (jTljl we obtain for the decay — > e~X\ the end points 
Ef a = 63 GeV, E™ x = 153 GeV, and for e L -> e~x\ the limits Ef° = 38 GeV, 
P™ ax = 93 GeV. In this scenario the SrCl background is sizeable, and negative beam 
polarisation might be needed in order to clearly determine the end points of the electron 
energy distributions. (Note that the large distortion in the energy spectra for tRth is 
a result of the reconstruction process, which is devised for the e^e-h and crCr signals.) 
The electron energy distributions in e^e^ production display the end points at the 
expected energies Ef n = 32 GeV, P x max = 166 GeV, Ef n = 12 GeV and P 2 max = 60 
GeV. For an integrated luminosity of 100 fb -1 , these end points may be hard to observe 
even with positive beam polarisation, because of the small cross section and the poor 
statistics for this process. However, for larger luminosities the end points could be 
observed more clearly (see Appendix iBl). 

5 Conclusions 

In this paper we have focused on the experimental test of the Majorana nature of the 
neutralinos, studying selectron pair production in e~e~ scattering. These processes 
can only take place if, as it is predicted in most SUSY models, the neutralinos are 
Majorana particles. Motivated by this issue, we have demonstrated that it is possible to 
reconstruct the selectron masses in the processes e~e~ — > e^e^, eRtR, using a relatively 
rare channel with one selectron decaying to e~x? an d the other one decaying via e — > 
e ~X2 ~^ e ~Xiff- The reconstruction can be done using as input the 4-momenta of the 
detected particles, the CM energy and the Xi an d X2 masses, which are assumed to be 
known. In addition, the reconstruction procedure allows for a determination of all final 
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Ej (GeV) 

Figure 13: Kinematical distribution of E 1 in scenario 2, for unpolarised beams. 
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Figure 14: Kinematical distribution of E 2 in scenario 2, for unpolarised beams. 

state momenta, which in turn has other potential applications [15J. We have illustrated 
our results in two mSUGRA scenarios. In the first scenario we have considered the 
channel // = /i + yU~, and in the second one the channel // = qq. In the calculations 
we have taken into account ISR, beamstrahlung and particle width effects, and we 
have performed a simple simulation of the energy resolution of the detector. We have 
included the main background from e^ei, production as well. We have shown that 
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for the kinematically allowed processes the masses can be successfully reconstructed in 
both scenarios, yielding a clear peak in the case of e^L production and a smaller peak 
for production. The examination of the electron energy spectra shows that the 

produced particles are scalars, making it obvious that the observed signal corresponds 
to selectron pair production. Additionally, the electron energy distributions display 
end points at the energies expected for e^e-h and crCr production in each of the cases 
analysed, providing independent determinations of the selectron masses and further 
supporting the evidence for selectron pair production. Obviously, the hypothesis of 
selectron pair production can also be confirmed by comparing the reconstructed masses 
with the selectron masses precisely measured in other processes [TfllfTT] . 

In e~e~ collisions, selectron pair production requires the exchange of a Majorana 
neutralino in the t channel, but does not prove that all the neutralinos are Majorana 
particles. However, if the neutralino mixing matrix is precisely reconstructed in other 
experiments pHHHUEH], the measurement of the e^L and tRtR cross sections may 
allow us to identify the contributions of the different neutralino mass eigenstates. In 
mSUGRA scenarios as those here considered, where the two lightest neutralino mass 
eigenstates are gaugino-like, crCr production is mainly mediated by the exchange of 
Xi, with a smaller contribution from the opposite occurring for e^-L production. 
Then, in these scenarios the measurement of one or both cross sections can prove 
the exchange of Majorana neutralinos Xi and X2 (the analogous can be done in more 
general scenarios where the reverse situation occurs, that is, when Xi is dominated by 
the wino component and X2 * s bino-like). The contribution to the cross sections of the 
neutralinos with large higgsino-like components (^3 and X4 m our case) is expected to 
be smaller than the uncertainty in the measurement of the cross sections. 

It is worth noting here that in addition to selectron pair production in e~e~ col- 
lisions, the processes e + e~ — > XiXi ~~* with i = 2,3,4 and I = e, /z, may 
also shed light on the Dirac or Majorana nature of the neutralinos [3 El El- The 
energy distributions of the final state charged leptons are sensitive to the Dirac or 
Majorana nature of the decaying neutralino Xi 5 allowing for the construction of lep- 
ton energy asymmetries, which vanish for a Majorana x® but may have a nonzero 
value if X4 is a Dirac particle. Unfortunately, these processes suffer from severe back- 
grounds from slepton pair production e + e~ — > 11 — > l + l~XiXi, as well as chargino pair 
e + e~ — > xtxj ~^ l + l 'vfiiXiXi and W + W~ production. Moreover, for the case of X1X3 
and XiXi the production cross sections and decay branching ratios are in general small. 
If these difficulties can be overcome, these processes may give further evidence of the 
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Majorana or Dirac nature of X3 an d xl- 

In this work we have taken as example the TESLA collider with a CM energy of 
500 GeV, but the analysis presented can be applied to the proposed TESLA upgrade 
to 800 GeV and other future colliders with larger CM energies, like CLIC. In this case, 
with the larger energies and luminosities the processes studied will be relevant in a 
larger region of parameter space, which is not kinematically accesible at 500 GeV. Due 
to the t-channel nature of these processes, their cross sections remain sizeable at TeV 
energies and hence allow us to investigate the Majorana nature of the neutralinos in 
the case of a heavier SUSY spectrum too. Finally, let us remark the importance of the 
reconstruction of the selectron masses on its own. This novel reconstruction technique 
can also be used, for instance, in selectron and smuon pair production in e + e~ collisions. 
It provides a new and independent measurement of the masses of these sfermions, and 
can be used as a tool for the analysis of the angular distributions in the production 
and decay of these sfermions, and for the study of their spin [11]. In addition, the 
mass reconstruction may be used in order to discriminate between these and other 
SUSY signals. For larger CM energies, the same method can be applied to squark pair 
production, with or without the help of flavour tagging. 
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A Lagrangian 
A.l Mass matrices 

At low energies, the charged slepton mass term in the weak eigenstate basis is 

£| ass = - l -E'^M\E' , (A.l) 
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with the most general charged slepton mass matrix 

/ M 2 - M? 

M l = I wf" Jo LR I • ( A -2) 




Omitting the flavour dependence, the above submatrices are given by jHEl 

M\ ll = ^ vl Y e *Y e T + Ml - ^ m 2 z cos2(3 (1-2 sin 2 9 W )1, 

M\ rr = l -vlYjY: + Ml R -m 2 zC os2P sin 2 ^ IL , 

M% LR = Ml L = ^=(Y-y-,^=Y : , (A.3) 

where Y e are the charged lepton Yukawa couplings, and M^ R the soft breaking 
scalar masses for the left and right handed sleptons, respectively, Y^ the soft trilinear 
terms and 1 is the 3x3 identity matrix in flavour space. The trilinear couplings can 
be decomposed as (Y A )ij = Y^, with no summation over The slepton mass 
matrix is diagonalised by a 6 x 6 rotation, 

(Mp) 2 = R,£ M| Fig t , (A.4) 

with E = R^E' the mass eigenstates. Neglecting the mixing between generations, 
the matrices in Eqs. (JA.2IA.3J) are diagonal in flavour space, and only LR mixing is 
present. Hence, the former relation between mass and weak interaction eigenstates can 
be simply expressed for each flavour as 

l = R l l\ (A.5) 

where R l now denotes a 2 x 2 matrix, with I — e, [A, r. Using the conventions of 
Ref. [201 5 an d assuming that the matrices in Eqs. (I A. 3 J) are real, 




sin Uj 



cos 9j sin 9 




(A.6) 



with m ( - < raj- . 



Squark-mediated interactions have a sub-dominant role in the parameter space ana- 
lysed. Similarly to the slepton case, the 6x6 squark mass matrices can be diagonalised 
by a unitary matrix Rq, 

(M^) 2 = Rq M% Rj, Q = U,D (A.7) 
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so that the relation between the weak (primed) and mass eigenstates is given by 



q = R d q' . 



(A.8) 



The four Majorana neutralinos x° are mixtures of weak interaction eigenstates 
(bino, neutral wino and neutral higgsinos). In the basis where (ip°) T = (B, W®, H®, H^) 7 , 
the mass term is 

£™o ass = ™ (^°) T Mvo ^° + H. c. , (A.9) 
2 

and the neutralino mass matrix can be written as 



mi — mz s'md\y cos/3 mz s'mOw sin/3 N ^ 

ra-i mz cos^vv cos/3 — mz cosBw sin/3 

-mz sin^vi/ cos/3 mz cos^w cos/3 — \i 

y mz sin#vy sin/3 — mz cos#w sin/3 — /i 

(A.10) 

where m 12 are the soft gaugino masses. This matrix is diagonalised by 



/ 



N* M 9 o N- 1 = M A T 9 . 



(A.ll) 



A. 2 Interaction terms 



In this section, we list the interactions 0EE1E3EE1 relevant for the calculation of the 
production and decay cross sections. The fermion-sfermion-neutralino terms are given 
by 



C 



km 



Si x 



~o 

3 



{C[)ijPL + (C R )ijP R 



f + fif 



(A.12) 



For the case of charged (s)leptons, the C couplings are 



V2 



(N* + tan 6 W N*) R l a - Y e N* R\ 2 , 



-gy/2twxB w N jX R l i2 - Y e N j3 R l a 



(A.13) 



with I = e, fi, t. In our computations we have considered the limit where Y e = 
diag(0, 0, h T ), so that the second term on the r.h.s. of the above equations is only 
present for the case of the taus. 

For the calculation of the decay x° — * Xi<?<? we include squark-mediated inter- 
actions, though they have a minor importance. In this case, the couplings have a 
slightly more cumbersome expression, which involves both quark and squark rotation 
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matrices. However, in this work, the effects of quark and squark flavour mixing are 
not relevant, and one can safely ignore them. Moreover, as in the lepton sector, the 
Yukawa couplings of the two first generations can be neglected. In this limit, where 
Y u = diag (0, 0, h t ), Yd = diag (0, 0, hb) and Vckm — 1, squark mixing is purely LR, 
and can be parametrised for every family as in Eqs. (|A.5| IA.6j) . With u — u, c, t and 
d = d, s, b, the C couplings read 

= -^(N; 2 + ^Un6 w N* 1 )Rf 1 -Y u N* 4 Rf 2 , 

iC%)ij = gV^^tanOwN^Rl-YuN^Rl, 
{C d L ) i3 = ^(N* 2 -~tan9 w N* 1 )Ri-Y d N* 3 Ri, 

(C d R )ij = -gV2^n9 w N 0l Ri-Y d N^Ri. (A.14) 

The Z-neutralino-neutralino interactions are parametrised by 



£z nn = Z, [xh» {DtP L + D%P R ) jg] , (A.15) 



where 



D% = \ (N U N* 4 - N i3 N* 3 ) , 

D% = -ipi)*. (A.16) 

Higgs-mediated interactions play a marginal role in the processes here analysed, and 
are only non-negligible for final states with // = bb, t + t~ . Assuming no CP violation 
in the Higgs sector (so that there is no mixture between the CP-even and the CP-odd 
states), the interaction of the neutralinos with the neutral CP-even physical Higgs can 
be written as 



£tf°x?x? = _ 7/| ( H ° cos a ~ h ° sin a ) X° [Qij p L + Q^Pr] X°j 

+ 4f (H° sin a + h° cos a) [S l3 P L + S&P R ] X° , (A.17) 



where a is the CP-even Higgs mixing angle, and as usual m h o < m H o. In Eq. (j A . 1 7|) 
we have the couplings 

% = ^JVS-tan^JVJa + ^j), 

S tj = liN^-tanOwN^N^ + ii^j). (A.18) 
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The interaction of the pseudoscalar with the neutralinos can be written as 



A f ° [ {Q tj sin (3 - Sij cos (3) P L - (Q* sin /? - 5* cos /3) P fl ] x 



.-,o 



(A.19) 

with the couplings as in Eq. (|A.18J) . The Higgs-fermion-fermion Lagrangian can be 
expressed, omitting the flavour dependence, as 

C Hff = -Yf ag [ (H° sin a + h° cos a) uu-i A cos (3 uj 5 u] 

_ydiag j ^0 cog a _ h gin a \ _ i A gin p J 7; . rf j 

-y e diag [ (H° cos a - h° sin a) M - i A sin /? l l5 l] . (A.20) 

The iff/ interactions in the above equation are proportional to the fermion Yukawa 
couplings, and hence are suppressed in all cases, the exception being the top and bottom 
quarks and the r. Finally, the Zff couplings are the standard ones, 



C 



zff 



9 



cos 9 



w 



l( - Q } sin 2 9 W )P L - Qt sin 2 6 W P R f , 



(A.21) 



where l( and Qf are, respectively, the weak isospin and charge of fermion /. 
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B Effect of beam polarisation 



The TESLA design offers the possibility of electron polarisation up to ±80%, which 
provides several advantages for our study, reducing the crc^ background and enhancing 
the or crCr signals. The cross sections for selectron pair production with lon- 
gitudinally polarised beams can be related to the unpolarised cross sections in a very 
simple way. Since selectron mixing is negligible, for initial electrons with definite helic- 
ity the only non-vanishing amplitudes for selectron pair production are e^e^ — > e^eL, 
e R e R ~ * an d ZrZl ~* &R&L- This allows to write the cross sections for arbitrary 
polarisations Pi, P 2 as 

(1-P 1 )(1-P 2 ) C7g igi (0,0), 
(1 + P0(1 + P 2 ) * SrSh (0,0) , 
(l-P 1 P 2 )a iRiL (0,0). (B.l) 

For the sake of simplicity, throughout this article we have plotted cross sections for 
unpolarised beams. For the case of polarised electrons, the cross sections are straight- 
forward to obtain, using Eqs. (|B.lj) . However, it is very illustrative to plot some of the 
distributions presented in Section 0] for the case of polarised beams, making apparent 
the enhancement of the signals and the reduction of CrCl production. We only consider 
the cases Pi = P 2 = —0.8 and Pi = P 2 = 0.8. The use of a polarisation Pi = — P 2 = 0.8 
does not offer any advantage for our study. 

B.l Negative beam polarisation 

Negative beam polarisation Pi = P 2 = —0.8 enhances the e^L signal by a factor of 
3.24, reducing the tRtR signal by a factor of 0.04 and the tRth background by 0.36. 
The effect on the mass distributions and electron energy spectra for scenario 1 can be 
seen in Figs. fTKHTTl For scenario 2, the corresponding cross sections are depicted in 
Figs. HHH13 



^e i e i (Pl,P 2 ) = 
&e R e R {Pl, P 2 ) = 
a e R e L {Pl-,P2) = 
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Figure 15: Reconstructed selectron mass in scenario 1, for P\ = P 2 = —0.8. 
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Figure 17: Kinematical distribution of E 2 in scenario 1, for Pi = P 2 = — 0. 




300 



m- (GeV) 



Figure 18: Reconstructed selectron masses in scenario 2, for Pi = P 2 = —0.8. 
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Figure 19: Kinematical distribution of E\ in scenario 2, for P\ = P 2 = —0.8. 




Figure 20: Kinematical distribution of E 2 in scenario 2, for P\ = P 2 = —0.8. 



31 



B.2 Positive beam polarisation 



Positive beam polarisation P\ = P 2 = 0.8 enhances the BrCr signal by a factor of 
3.24, reducing the e^L signal by a factor of 0.04 and the e^e^ background by 0.36. 
In scenario 1, positive beam polarisation does not offer any advantage, and we only 
present the plots for scenario 2, in Figs. l2lT|23l 
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Figure 21: Reconstructed selectron masses in scenario 2, for P\= P 2 = 0.8. 
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Figure 22: Kinematical distribution of E x in scenario 2, for P\ = P 2 = 0. 
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Figure 23: Kinematical distribution of E 2 in scenario 2, for Pi = P 2 = 0.8. 
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